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Abstract 

The problem of non-perturbative description of stationary flames with arbi- 
trary gas expansion is considered. A general method for deriving equations 
for the flame front position is developed. On the basis of the Thomson circu- 
lation theorem an implicit integral of the flow equations is constructed. With 
the help of this integral the flow structure near the flame front is determined, 
and a closed system of equations for the flame front position and a set of 
auxiliary quantities is obtained. This system is shown to be quasi-local, and 
its transverse representation is found. This representation allows reduction of 
the system to a single equation for the flame front position. The developed 
approach is applied to the case of zero-thickness flames. 
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I. INTRODUCTION 



The process of flame propagation presents an extremely complicated mathematical prob- 
lem. The governing equations include the nonlinear flow equations for the fuel and the 
products of combustion, as well as the transport equations governing the heat conduction 
and species diffusion inside the flame front. Fortunately, in practice, an inner flame scale 
defined by the latter processes is large compared to the flame front thickness, implying that 
the flame can be considered as a gasdynamic discontinuity. The initial problem is reduced 
thereby to a purely hydrodynamic problem of determining the propagation of a surface of 
discontinuity in an incompressible fluid, the laws of this propagation being determined by 
the usual Navier-Stokes and continuity equations complemented by the jump conditions at 
the surface, expressing the mass and momentum conservation across the flame front. The 
asymptotic methods developed in Refs. allow one to express these conditions in the 

form of a power series with respect to the small flame front thickness. 

Despite this considerable progress, however, a closed theoretical description of the flame 
propagation is still lacking. What is meant by the term "closed description" here is the 
description of flame dynamics as dynamics of the flow variables on the flame front surface. 
Reduction of the system of bulk equations and jump conditions, mentioned above, to this 
"surface dynamics" implies solving the flow equations for the fuel and the combustion prod- 
ucts, satisfying given boundary conditions and the jump conditions at the flame front, and 
has only been carried out asymptotically for the case — > 1, where 9 is the gas expansion 
coefficient defined as the ratio of the fuel density and the density of burnt matter Ijl],^-^. 

Difficulties encountered in trying to obtain a closed description of flames are conditioned 
by the following two crucial aspects: 

(1) Characterized by the flow velocities which are typically well below the speed of sound, 
deflagration represents an essentially nonlocal process, in the sense that the flame-induced 
gas flows, both up- and downstream, strongly affect the flame front structure itself. A seeding 
role in this reciprocal process is played by the Landau-Darrieus (LD) instability of zero 
thickness flames A very important factor of non- locality of the flame propagation is the 
vorticity production in the flame, which highly complicates the flow structure downstream. 
In particular, the local relation between pressure and velocity fields upstream, expressed by 
the Bernoulli equation, no longer holds for the flow variables downstream. 

(2) Deflagration is a highly nonlinear process which requires an adequate non-perturba- 
tive description of flames with arbitrary values of the flame front slope. As a result of develop- 
ment of the LD-instability, exponentially growing perturbations with arbitrary wavelengths 
make any initially smooth flame front configuration corrugated. Characteristic size of the 
resulting "cellular" structure is of the order of the cutoff wavelength Ac ~ 20Lf given by 
the linear theory of the LD-instability 0; Lf is the flame front thickness. The exponen- 
tial growth of unstable modes is ultimately suppressed by the nonlinear effects. Since for 
arbitrary 9 the governing equations do not contain small parameters, it is clear that the 
LD-instability can only be suppressed by the nonlinear effects if the latter are not small, 
and therefore so is the flame front slope. 

The stabilizing role of the nonlinear effects is best illustrated in the case of stationary 
flame propagation. Numerical experiments on 2D flames with ^ = 6 — 8 show that even 
in very narrow tubes (tube width of the order Ac), typical values for the flame front slope 
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are about 1.5 — 2.0 0. Nonlinearity can be considered small only in the case of small gas 
expansion, 6-^1, where one has the 0(^ — 1) estimate for the slope, so that it is possible 
to derive an equation for the flame front position in the framework of the perturbation 
expansion in powers of {9 — 1) [|-|^. 

This perturbative method gives results in a reasonable agreement with the experiment 
only for flames with 6 < 3, propagating in very narrow tubes (tube width of the order Ac), 
so that the front slope does not exceed unity. Flames of practical importance, however, 
have 6 up to 10, and propagate in tubes much wider than Ac. As a result of development 
of the LD-instability, such flames turn out to be highly curved, which leads to a noticeable 
increase of the flame velocity. In connection with this, a natural question arises whether 
it is possible to develop a non-perturbative approach to the flame dynamics, closed in the 
sense mentioned above, which would be applicable to flames with arbitrary gas expansion. 

A deeper root of this problem is the following dilemma: On the one hand, flame propaga- 
tion is an essentially non-local process [see the point (1) above], on the other, this non-locality 
itself is determined by the flame front configuration and the structure of gas flows near the 
front, so the question is whether an explicit bulk structure of the flow downstream is nec- 
essary in deriving an equation for the flame front position. In other words, we look for an 
approach which would provide the closed description of flames more directly, without the 
need to solve the flow equations explicitly. 

The purpose of this paper is to develop such approach in the stationary case. 

The paper is organized as follows. The flow equations and related results needed in our 
investigation are displayed in Sec. [II A| . A formal integral of the flow equations is obtained 
in Sec. [II B| on the basis of the Thomson circulation theorem. With the help of this integral, 
an exact flow structure near the flame front is determined in Sec. pT[ Using the jump 
conditions for the velocity components across the flame front, this result allows one to write 
down a system of equations for the flow variables and the flame front position, which is 
closed in the above-mentioned sense. This is done in the most general form in Sec. [IV A . 
It is remarkable that this system turns out to be quasi-local, which implies existence of the 
transverse representation derived in Sec. IV B[ . The developed approach is applied to the 
particular case of zero-thickness flames in Sec. |IV C| , where a single equation for the flame 
front position is derived. The results obtained are discussed in Sec. 0. 



II. INTEGRAL REPRESENTATION OF FLOW DYNAMICS 

As was mentioned in the point (1) of Introduction, an important factor of the flow 
non-locality downstream is the vorticity production in curved flames, which highly compli- 
cates relations between the flow variables. In the presence of vorticity, pressure is expressed 
through the velocity fleld by an integral relation, its kernel being the Green function of the 
Laplace operator. It should be noted, however, that the jump condition for the pressure 
across the flame front only serves as the boundary condition for determining an appropriate 
Green function, being useless in other respects. Thus, it is convenient to exclude pres- 
sure from our consideration from the very beginning. The basis for this is provided by the 
well-known Thomson circulation theorem. Thus, we begin in Sec. |II A[ with the standard for- 



mulation of the problem of flame propagation, and then construct a formal implicit solution 
of the flow equations with the help of this theorem in Sec. [ilB. 
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A. Flow equations 



Let us consider a stationary flame propagating in an initially uniform premixed ideal fluid 
in a tube of arbitrary width b. To make the mathematics of our approach more transparent, 
we will be dealing below with the 2D case. Let the Cartesian coordinates (x, z) be chosen so 
that 2;-axis is parallel to the tube walls, z = — oo being in the fresh fuel. It will be convenient 
to introduce the following dimensionless variables 

iv, = i^/R, z/R) , {w, u) = {yjUi , yjUi) , 
n= (P-Pf)/p_f/f2, 

where Uf is the velocity of a plane flame front, Pf is the initial pressure in the fuel far 
ahead of the flame front, and R is some characteristic length of the problem (e.g., the cut- 
off wavelength). The fluid density will be normalized on the fuel density p_. As always, 
we assume that the process of flame propagation is nearly isobaric. Then the velocity and 
pressure flelds obey the following equations in the bulk 

I-' « 

where (C15C2) = iv^O^ {vi,V2) = {w,u), and summation over repeated indices is implied. 

Acting on Eq. (H) by the operator Skid/dC^u where Sik = —£ki, £12 = +1, and using 
Eq. (|1]), one obtains a 2D version of the Thomson circulation theorem 

..,| = 0. (3) 

where 

_ du dw 
drj 

According to Eq. (|]), the voticity a is conserved along the stream lines. As a simple 
consequence of this theorem, one can flnd the general solution of the flow equations upstream. 
Namely, since the flow is potential at ^ = —00 {u = V = const, w = 0, where V is the velocity 
of the flame in the rest frame of reference of the fuel), it is potential for every ^ < f{ri), f{rj) 
denoting the flame front position. Therefore, 



hoo 

u = 

n=— 00 



Y: w(")exp -(|n|e + mr/) , (4) 
i=— 00 ^ 

w = H{u- u(°)) , (5) 
where the linear Hilbert operator H is deflned by 

H explikrj) = ix{k) exp^ikrj) , 7^ , (6) 

and 

, . / +1, X > 0, 
^(^)= -1, x<0. 
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The boundary condition w = at the tube walls (r/ = 0, b) implies that the coefficients u^"'^ 
are all real, 

= . (7) 

It will be shown in the next section how the Thomson theorem can be used to obtain a 
formal integral of the flow equations downstream. 



B. Integration of the flow equations 

Consider the quantity 

/In T 

s 

where r is the distance from an infinitesimal fluid element ds to the point of observation (, 
= {Q — (i)"^, and integration is carried over S = {fj,^ : C, > f{v)}- Taking into account 
Eq. (|l]), one has for the divergence of Oj : 

wrl^' (^) (^) (^'^) = ~ / ^'^^^^ ' 

where di = d/dQ, A boundary of S, and dli its element. 

Next, let us calculate SikOkSimdiam- Using Eq. (|^), we find 

„ „ d ( da2 daA d ( dai f\nr , \ d'^ai 

£lkOk£lmOiam = TTT KT =7^ — 7: / -^Vidk 



\ dr] J drj \ drj •/ 27r j d^"^ 

d /■ Inr 



A 



Aai - T- f ^Vidli . 



drj 

Analogously, 

£2kdkeimdiam = - Aa2 - J "^^^'^^^ • 

^ A 

Together, these two equations can be written as 

/In r 
-^Vkdlk ■ 

A 

Substituting the definition of into the latter equation, and integrating by parts gives 

Inr 



-EikOkEi^diA -di f -^Vkdlk 

J 271 

A 

f Inr f Inr f Inr . . 

SikOk J —eimVmdh - Oi j -^Vkdlk - SikOk J i^^ds . (9) 



A A E 
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The first two terms on the right of Eq. @ represent the potential component of the fluid 
velocity, while the third corresponds to the vortex component. The aim of the subsequent 
calculation is to transform the latter to an integral over the flame front surface. To this end, 
we will decompose S into elementary ds as follows. 

Let us take a couple of stream lines crossing the flame front at points {i],f{ri)) and 
[rj + Ari,f{ri + Arj)) (see Fig. p. Consider the gas elements moving between these lines, 
which cross the front between the time instants t = and t = At. During this time interval, 
these elements fill a space element Sq adjacent to the flame front. For sufficiently small 
Ar], At, the volume of Eg 



As 

where 



Ar] f'Ar] 

w+At u.At 



- f'w+)Ar]At = vlNAT]At 



the subscript "+" means that the corresponding quantity is evaluated just behind the flame 
front, i.e., for ^ = /(r/) + 0, and f" = fj+rij is the normal velocity of the burnt gas, rii being 
the unit vector normal to the flame front (pointing to the burnt matter). After another time 
interval of the same duration At, the elements move to a space element Si adjacent to Sq. 
Since the flow is incompressible. Si is of the same volume as Sq. Continuing this, the space 
between the two stream lines turns out to be divided into an infinite sequence of S's of the 
same volume, adjacent to each other. Thus, summing over all Arj, the third term in Eq. (^ 
can be written as 

-^9, / dlvla+K{7],^,f]), (10) 



F 



2 

where F denotes the flame front surface (the front line in our 2D case), 

1 °° 1/2 

K{7], r/) = - lim E 1^ {(r/ - H(r/, riAt))^ + (^ - S(r/, riAt))^} At 

+ 00 

= -/ dtln{{v-mt)f + {^~E{v,t)fy^\ (11) 



and (H(?7, t), S(f/, t)) trajectory of a particle crossing the point {fj, fifj)) at t = 0. 
Substituting expression (|^) into Eq. @) gives 

Vi = Eikdk j ^^eimVmdli -di J ^-^Vkdk - J dl vla+K{r],^,f]) . (12) 

A A F 

This representation of the flow velocity downstream will be used in the next section to 
determine the structure of the vortex mode near the flame front. 
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III. STRUCTURE OF THE VORTEX MODE 



To determine the flame front dynamics, it is sufficient to know tlie flow structure near 
the flame front. As to the flow upstream, it is described by Eqs. (^, (|) for all ^ < f{rj). 
Given the solution upstream, velocity components of the burnt gas at the flame front can be 
found from the jump conditions which express the mass and momentum conservation across 
the front. On the other hand, these components are required to be the boundary values 
(for ^ = + 0) of the velocity field satisfying the flow equations. As was shown in the 
preceding section, the latter can be represented in the integral form, Eq. ([T^) . Any velocity 
field can be arbitrarily decomposed into a potential and vortex modes. The former has the 
form analogous to Eqs. (11), (|), namely 



u 



^= Y: M(")exp]^(-|n|e + ^n77)f , (13) 



TT 



= -H{uP - m(°)) . (14) 
Our strategy thus will be to use the integral representation to determine the near-the-front 



structure of the vortex mode described by the last term in Eq. ([121) . 

Equation (|T2|) reveals the following important fact. Up to a potential, the value of the 
vortex mode at a given point (?], of ^"^^ downstream is determined by only one point 
in the range of integration over (?7,t), namely that satisfying 

H(77,t) = r7, 2(77, t)=e. (15) 

This is, of course, a simple consequence of the Thomson theorem underlying the above 
derivation of Eq. ([T^). It can be verified directly by calculating the rotor of the right hand 
side of Eq. ([121) : contracting this equation with Eimdm, using eik£im = Skm, and taking into 
account relation 

Alnr = 2tt6^^\C - , 



one finds 

F OF 



le^md^e^kdk ! dlvla+{fi)K{r],i,f^)= I dt I dl vla4fi)6{r] -R{fi,t))6{^ - E{fi,t)) . (16) 



Since 5{x) = for a; 7^ 0, the product of 5-functions picks the point (|1^) out of the whole 
range of integration in the right hand side of Eq. ([l^). 

Now, let us take the observation point (^,?7) sufficiently close to the flame front, i.e., 
^ ^ f{ri), > f{r])]. In view of what has just been said, the vortex component for such 
points is determined by a contribution coming from the integration over fj, t near the flame 
front, which corresponds to small values of t. Integration over all other f/,t gives rise to a 
potential contribution. 

The small t contribution to the integral kernel K[ri,C,,'fi) can be calculated exactly. For 
such t's, one can write 

R{fl,t) ^fl + w+{fi)t, E{fi,t) ^ f{fi)+u+{fi)t. (17) 
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Let the equality of two fields fiirj,^), ^2{j]iC) up to a potential field be denoted as ipi = ip2- 
Then, substituting Eq. ( |T7|) into Eq. (|l3), and integrating gives 



+00 


1 1/2 
= - J dt\n {vlf - 2rv+t + r^}^ ^ 



1/2 







- J dy In \ - 

rvi V 



1/2 



\ 



rv I 



arctan 



-y + y\n{y'^ - 



1/2- 



(rv+)^ 



v+to- 



(18) 



Here v+ denotes the absolute value of the velocity field at the fiame front, and to is assumed 
small enough to justify the approximate equations (|T7|) . 

As we know, the only point in the range of integration over fj, t that contributes to the 
vortex mode is the one satisfying Eq. (pIS]) or, after integrating over t, 



[v - - - fifi)]w+ifi) = . 



(19) 



The distance r between this point and the point of observation tends to zero as the latter 
approaches to the fiame front surface. Thus, taking — /(//)] small enough, one can make 
the ratio to/r as large as desired; therefore, the right hand side of Eq. (KW) is = 



rv I 



rv4 



Tc / rv+ 

— h arcsm 

2 y rv+ 



rv^ rv+ 

rt>+ rt>_|_ 



Inr ^ +TIC, 



(20) 



where "TIC" stands for "Terms Independent of the Coordinates" {i],C,)- Denoting 

rv_|_ 



n 



we finally obtain the following expression for the integral kernel 

K{7],i,fj) = |Vl -1^2 (^^ + aTcsmi?j +filn^| + TIC. 

In order to find the vortex mode of the velocity according to Eq. (|I2|) we need to calculate 
derivatives of K. Using relation 
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one easily obtains 

dK . I (vi 7r/2 + arcsinil Vi+ 

— — = \l , Inr. 22 

Equation (p2D can be highly simplified. Consider the quantity 

d { ( Ti ^Vi> \ arcsinf] — 7r/2 t>i . , 1 
T = —{\--VL—\ ' +—\nr\. (23) 



First, we calculate 



d (vi QYi±\ = -r-^ -—— = — (24) 



Second, we note that the vector 



satisfies 

= 1, /?ifj+ = 0, 

i.e., Pi is the unit vector orthogonal to v+ . In addition to that, (3i changes its sign at the 
point defined by Eq. (plQl). Therefore, the derivative of jSi, entering T, contains a term with 
the Dirac (5-function. However, this term is multiplied by (arcsini7 — vr/2), which turns into 
zero together with the argument of the 5-function. Thus, using Eqs. ([2l|) , (^) , (^) , one finds 



fl^ arcsin Q — 7t/2 ( Vi ^ Vi< ^ 

T = ' + - - ^— 

r V 1 - 1]2 \r v+ ^ 

We conclude that the term 



1 ^ arcsin Q — 11/2 



on n 

7^ + - = 0- 



f Ti arcsin 11 - 7r/2 Vi+ 

\l — , mr 

\r v+J v^r^TF v+ 

in the integral kernel corresponds to a pure potential. Therefore, we can rewrite Eq. (E2|) as 



^ ^ h r^^^A / - (26) 

Finally, substituting this result into Eq. ([T^, noting that the vector SkiPk is the unit vector 
parallel to f j+ if eikriVk+ > 0, and antiparallel in the opposite case, we obtain the following 
expression for the vortex component, v^, of the gas velocity downstream near the fiame front 
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Vi = J dl xi^pqrpVq+ f^^^^'^ . (27) 

F ^ 

Having written the exact equality in Eq. (P7| ) we take this equation as the definition of the 
vortex mode. As a useful check, it is verified in the appendix that the obtained expression 
for satisfies 



IV. CLOSED DESCRIPTION OF STATIONARY FLAMES 

After we have determined the near-the-front structure of the vortex component of the 
gas velocity downstream, we can write down a closed system of equations governing the 
stationary flame propagation. As was explained in the Introduction, the term "closed" 
means that these equations relate only quantities defined on the flame front surface, without 
any reference to the flow dynamics in the bulk. This system consists of the jump conditions 
for the velocity components at the front, and the so-called evolution equation that gives 
the local normal fuel velocity at the front as a function of the front curvature. These 
equations (except for the evolution equation) are consequences of the mass and momentum 
conservation across the flame front. In Sec. [IV A| , we obtain the closed system in the most 



general form, without specifying the form of the jump conditions, and then apply it to the 
case of zero thickness flames in Sec. IV C. 



A. General formulation 

First of all, we need to find the "on-shell" expression for the vortex component, i.e., its 
limiting form for ^ — > /(r/) + 0. To this end we note that in this limit, x{^ik^i'^k+) xiv^v)^ 
therefore, Eq. (|27|) gives 

v^^ = ldlxiv-vf-^^- (28) 

F + 

The total velocity field downstream is the sum of the potential mode having the form 
Eqs. (p!3|) , ([I^) , and the vortex mode. Let us denote the jump [vi+{ri) — Vi-^rj)] of the gas 
velocity Vi across the flame front as [vi]. Here v±i{r]) = Vi{ri, fij]) ± 0). Then we can write 

The jumps [vi\ (as well as cr+) are quasi-local functions of the fuel velocity at the flame 
front, and of the flame front shape. Two equations (p9|), together with Eqs. (H), (0), and 
the evolution equation, = v^{f), form a closed system of five equations for the five 
functions Vi-{ri),vf_^{ri), and f{ri). 

It should be emphasized that this formulation implies that the potentials Vi and vf are 
explicitly expressed in the form of the Fourier series, Eqs. (^) and (|13[), respectively. Indeed, 
relations (^, (|I4D between the flow variables hold before the on-shell transition ^ — > /(//) 
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is performed, while Eq. (^) is formulated in terms of the on-shell variables v±i. Thus, the 
above system is in fact a system of equations for the front position f{ri) and two infinite sets 
of the Fourier coefficients u^'^\u^^\ 

However, the form of the integral kernel in Eq. ( ^9]) makes it possible to avoid this consid- 
erable complication, and to derive a much simpler formulation. This will be demonstrated 
in the next section. 



B. Transverse representation 

The aforesaid simplification is based on the fact that Eq. (|29D is quasi-local. To show 
this, we simply increase its differential order by one. Namely, differentiating this equation 
with respect to 77, taking into account relation 

= 2(5(x) , (30) 
and performing the trivial integration over f/, we get 



Given a point at the flame front, Eq. (|3TD relates components of the fuel velocity at this 
point and their derivatives at the same point. This property of quasi-locality implies ex- 
istence of the transverse representation of Eq. (|3T|). We say that a system of equations 
is in the transverse representation if all operations involved in this system (differentiation, 
integration) are only performed with the flrst argument (r/) of the flow variables. In other 
words, the ^-dependence of the flow variables in such system is purely parametric. 

Now, let us show how the system of Eqs. (||), (|1^, and ( pi]) can be brought to the 



transverse form. All we have to do is to express the full //-derivatives in terms of the partial 
ones. This is easily done using the continuity equation (|ip and the potentiality of the flelds 
ff as follows: 

(^)_^^'(^)_^ (^)_"^'(^)_ ^^^^ 

(34) 
(35) 

As to Eqs. (H), (fl4D , they are already transverse. Finally, the evolution equation has general 
form 

where F is a quasi-local function of its arguments, proportional to the flame front thickness, 
and therefore can also be rendered transverse. 



drj 

du\ 
drj 

dw\ 
drj 
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Thus, the complete system of governing equations can be written as follows: 



~a / ^ + M = ^ / ^ ^ 



df] drj drj drj + [u^]* 

dw ^,du , dwP ^,duP , (1 + F* + + H*) 



drj drj drj drj f ^ + [v'^Y 



w 



(*) 



y u-f'w = N + NF\ J 



The superscript "t" in these equations means that the corresponding quantity is to be 
expressed in the transverse form using Eqs. (|32D - (^), yet without setting ^ = in their 
arguments. The latter operation is displayed out the large brackets in (*). 

The meaning of the transformations performed is that the ^-dependence of the flow vari- 
ables is now irrelevant for the purpose of derivation of equation for the flame front position. 
Indeed, since all operations in the set of equations (*) are carried through in terms of rj 
only, one can solve these equations with respect to f{ri) under the large brackets signifying 
the on-shell transition. Furthermore, since the function f{ri) is itself .^-independent, we can 
consider all the flow variables involved in (*) ^-independent (since the resulting equation for 
/(r/) does not contain these variables anyway), and to omit the large brackets.^] 

C. Zero-thickness flames 

As an illustration, the developed approach will be applied in this section to the case of 
zero-thickness flames. 

In this case, the jump conditions for the velocity components have the form 

[«] = ^ , (36) 
H = -/'^ , (37) 



while the evolution equation 



1 . (38) 



We see that the jumps are velocity- independent, and F = 0. Also, it follows from these 
equations that 

[v^^] = 9-1, [v^] = 9^-1. 
Next, substituting the Fourier decomposition, Eq. (^, into evolution equation written as 

u- - f'w- = N , (39) 



^This reasoning has been repeatedly used in Refs. |]5|,|6|,|To| in deriving perturbative equations for 
the flame front position in the stationary as well as non-stationary cases. 
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and taking into account Eq. we get 
Integrating this equation over interval (0, b) gives 

b 

hu^^^ = [ dr]N = bV, or = V . 



It remains only to calculate the value of the vorticity at the flame front, as a function 
of the fuel velocity. This can be done directly using the flow equations (|ID,(0), see Ref. 0. 
With the help of Eqs. (5.32) and (6.15) of Ref. 0, the jump of the vorticity across the flame 
front can be written, in the 2D stationary 



cr 



9N \ ' ' ' N 

where 



Dw- + f'Du- + ^Df ] , (40) 



Differentiating the evolution equation (|39|), and writing Eq. (|40|) longhand, expression in the 
brackets can be considerably simplified 

1 (f'w y if'Yu' N' 

Dw^ + fDu_ + j^Df ^ w'_w_ + ^-l^ + + fu'_w_ + ^ 

= ^ + ^ ^ ' + + {u^-N) + —= uiu^ + w'_w. . (42) 

Since the flow is potential upstream, we obtain the following expression for the vorticity just 
behind the flame front 

g 

^+ = — opT^^'^^- + ^-^-) • (43) 

Using Eqs. (|3^) , (|33D in Eq. (|^), substituting the result into (*), and omitting the large 
brackets we arrive at the following fundamental system of equations 

^'_/V+(e-l) i^l^' = yP' - f'^vP 

-{9 - 1) {v iv' - /V) + c (c' + /V)} , (44) 

oo' + fv'-ie-l) (^^^ =^vP + fvP' 

-{0 - 1) {V {v' - /V) + fv')} "r+^JV/^^-T ' (45) 

u = H{v-V), (46) 
v-fuj = N, (47) 
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where v,uj,v^,uj^ are the ^-independent counterparts of the flow variables u,w,uP,w^, re- 
spectively, and $ = —Hd/drj is the Landau-Darrieus operator. 

It is not difficult to reduce the above system of equations to a single equation for the 
function fijf). For this purpose it is convenient to denote v = v — V . Then Eqs. (^61), (^T]) 
can be solved with respect to v to give 

v = {l- fH)-\N -V) = {1 + f'H + f'Hf'H + ...){N-V). (48) 

The two remaining Eqs. (H3|) and (ESI) are linear with respect to v^, . To exclude , we 



will use a complex representation of these equations. Namely, we multiply Eq. (^) by i and 
add it to Eq. (||): 

{v' + zuj')ii + zf) + {e-i) (^r-) = + ^^vm^ + 

-iO - 1) . . - W)} -"-rlllj^^^Zf^/^ . (49) 

It follows from the definition of the Hilbert operator, Eq. (^, that 

H"^ = -1. 

Therefore, 

(1 + iH){vP' + i^v^) = (1 + iH){l - iH)vP' = (1 - i^H^)vP' = . 



Thus, dividing Eq. (]49|) by (1 + if), acting by the operator (1 + iH) from the left, and 
taking into account that 

(1 + iH){v' + iuj') = (1 + iH){l + iH)v' = (1 + 2iH + i^H'^)v' = 2(1 + iH)v' , 

we obtain an equation for the fiame front position 



where {; is given by Eq. (^81). 

Equation (|50|) provides the closed description of the stationary zero-thickness fiames in 
the most convenient form. If desired, one can bring it to the explicitly real form by extracting 
the real or imaginary part.0 Account of the effects inside the thin flame front changes the 
right hand side of this equation to 0{e), e = L^/ R being the relative flame front thickness. 



^Both ways, of course, lead to equivalent equations: Acting by the operator H on the real part 
gives the imaginary one, and vice versa. 
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V. DISCUSSION AND CONCLUSIONS 



The results of Sec. |V| solve the problem of closed description of stationary flames. The 
set of equations (*) obtained in Sec. [IV B| gives a general recipe for deriving equations govern- 
ing the "surface dynamics" of the fields u,w,u^,w^, or, more precisely, their .^-independent 
counterparts, v,u!,v^,u!P, and the function f{r]) - the flame front position. Following this 
way, we derived an equation for the front position of zero-thickness flames, Eq. This 
equation is universal in that any surface of discontinuity, propagating in an ideal incom- 
pressible fluid, is described by Eq. (|50D whatever internal structure of this "discontinuity" 
be. The latter shows itself in the 0(£:)-corrections to this equation, where e is the relative 
thickness of discontinuity. In the case of premixed flame propagation, these corrections can 
be found using the results of Refs. following the general prescriptions of Sec. [IV B| . It is 



clear that independently of the specific form of e-corrections, the set (*) ultimately reduces 
to a single equation for the function f{ri), since equations of this set are linear with respect 
to the field . The latter can be eliminated, therefore, in the same way as in Sec. [IV C. 



It is interesting to trace the influence of the boundary conditions on the form of Eq. (pO|). 
By itself, expression ( [27[ ) for the vortex mode is completely "covariant", i.e., it has one and 
the same vector form whatever boundary conditions. The jump conditions Eqs. (0), ([37[) 
also can be rewritten in an explicitly covariant form as the conditions on the normal and 
tangential components of the velocity. It is the structure of the potential mode upstream and 
downstream, given by Eqs. (|p, (JTBD, respectively, which is directly affected by the boundary 
conditions for the velocity field. Thus, it is the boundary conditions which dictate the way 
the jump conditions appear in the first two equations of the set (*), as well as the form of 
the third and forth equations. 

Finally, it is worth to compare the results obtained in this paper and in Refs. where 
an equation describing stationary flames with arbitrary gas expansion was derived under 
assumption that there exists a local relation between the pressure field and a potential mode 
of the velocity downstream, expressed by the Bernoulli-type equation.^ This assumption 
was proved in Refs. 0,^] in the framework of perturbative expansion with respect to — 1) 
up to the sixth order. Comparison of Eq. ( [5(]| ) and Eq. (40) of Ref. now shows that 
this assumption is generally wrong (that the two equations are not equivalent can be easily 
verified, for instance, considering the large flame velocity limit investigated in detail in Sec. V 



of Ref. Ig]). As we saw in Sec. [II B| , the use of the Thomson theorem makes investigation 
of the pressure-velocity relation irrelevant to the purpose of deriving equation for the flame 
front position. 

The results presented in this paper resolve the dilemma stated in the Introduction in the 
case of stationary flames. There remains the question of principle whether it can be resolved 
in the general non- stationary case. 



^This assumption goes back to the work where it was introduced in investigating the nonlinear 
development of the LD-instabiUty. 
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APPENDIX: CONSISTENCY CHECK FOR EQ. (27) 



After a lengthy calculation in Sec. |T|, we obtained the following simple expression for 
the vorticity mode near the flame front 

< = J dl Xie^m^-^^ . (Al) 

As this important formula plays the central role in our investigation, a simple consistency 
check will be performed here, namely, we will verify that given by Eq. (^11) satisfies 

{duyOri - dwydO+ = (T+ . (A2) 

Contracting Eq. (|A1|) with Ekidk, and using relation (^), one finds 

f (TV' f 

SkidkVi = dl 6{epqrpVq+)ekiekmVm+ 2 = dl 6{epgrpVq+)vla+ . (A3) 



The argument of the 5-function turns into zero when the vectors and are parallel. 
Near this point, one can write 

£pqrpVg+ ^ rv+cp , 

where is the angle between the two vectors. On the other hand, the line element, dl, near 
the same point can be written as 

dl = d(j) = rd(f)— , 

sin ip 

simple geometric consideration shows, see Fig. 
Substituting these expressions into Eq. (|A3D, and taking into account relation 

6(ax) = - — :S(x) , 
\a\ 

we finally arrive at the desired identity 

d(j) 6{rv+(j))vla+ = / d(f) 6{(j))a+ = a+ . (A4) 



It should be noted in this respect that the identity Eq. ( |A2| ) is only a necessary condition 
imposed on the field v^. Playing the role of a "boundary condition" for the vortex mode, 
Eq. ( |A2| ) is satisfied by infinitely many essentially different fields, i.e., fields which are not 
equal up to a potential. It is not difficult to verify, for instance, that the velocity field defined 
by 

Vi = dl xiepgrpUq)—- 



also satisfies Eq. ( |A2| ) , and the difference v'^ — is essentially non-zero. 

By the construction of Sec. |T|, given by Eq. ( |A1|) is essentially the only field that 
satisfies the flow equations (|l|)-(@). 
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FIGURES 





FIG. 2. Near-the-front structure of the flow downstream, v" = nu" is the normal component 
of the velocity. Since the observation point {r],^) is close to the flame front, the stream line and 
the part of the front near this point can be considered straight. 
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